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Abstract 

Visualization of scalar and vector fields in the immersive environment (CAVE - Cave 
Automated Virtual Environment) is important for its application to radiation shielding 
researeh at NASA Eangley Research Center. A eomplete methodology and the 
underlying software for this purpose have been developed. The developed software has 
been put to use for the visualization of the earth’s magnetic field, and in partieular for the 
study of the South Atlantic Anomaly. The methodology has also been put to use for the 
visualization of geomagnetically trapped protons and eleetrons within Earth’s 
magnetosphere. 



Introduction 


Visualization of scalar and vector fields in the immersive environment (CAVE - Cave 
Automated Virtual Environment) is important for its application to radiation shielding 
researeh at NASA Eangley Research Center. Research has been performed for the 
development of a complete methodology and the underlying software for this purpose. 
Using the technique developed, a vector field ean be visualized in two modes: 

In mode one, the direetional eurved lines can be visualized (Eig.l.) in the surrounding 3- 
D space for any source(s) of a vector field. The size of the surrounding spherieal space 
for visualization of the direetional eurved line flux as well as the seale for this 
visualization ean be continually varied using the CAVE input eontroller . 

In mode two, the vector field at regularly spaeed points on the surfaee of a surrounding 
sphere for sources of the veetor field is represented correctly through the length and 
direction of the depicted veetors at these points (Fig.2.). As in mode one, the radius of the 
surface of the visualization sphere, and the seale for that visualization can be eontinually 
varied through the use of selected buttons on the eontroller . 

Customized modifications to the developed software have been introdueed for application 
to the visualization of the earth’s magnetic field, and in particular, for the study of the 
South Atlantie Anomaly (Fig.3.). 
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Fig. 1 : Field lines in the surrounding 3-D spherical region for 
a group of three electrical charges . The sizes, colors, and 
centers of the three spheres in the central region represent 
the magnitudes, polarity, and locations of the charges 
involved, respectively. 
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Fig. 2: Electrical field vectors on the surface of a surrounding 
spherical surface for the same three-charge group as 
described in the caption of Fig. 1 . 
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Fig. 3: Magnetic field of the earth in the surrounding 3-D 
spherical region. 
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Fig. 4: Geomagnetically trapped protons (red cloud) and electrons (white 
cloud) together with the magnetic field of the earth in the surrounding 3-D 
spherical region. 


A scalar field by itself can be visualized through the depiction of points with their 
number density proportional to the strength of the scalar field distribution in 3-D space. 
One other way the field can be visualized is through the depiction of small spheres in 3-D 
space with sizes proportional to the magnitude of the scalar field. The differentiation 
between positive and negative values of the field can be made through the use of two 
colors for plotted points or spheres, one color for the positive values of the scalar field, 
and the other color for the negative values. As an example of scalar field visualization, in 
Fig. 4, we show the visualization of the geomagnetically trapped protons and electrons 
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during solar maximum and minimum respectively, together with the earth’s magnetic 
field in the surrounding 3-D space of the earth. 

Project Description and Formulation 

Visualization of a Vector Field 


As stated in the introduction, a vector field can be visualized in two modes. In mode one, 
the directional curved lines are visualized in surrounding 3-D spherical space of chosen 
radius. To create, an approximately uniform density of these lines in the 3-D space, we do 
the following: 

The 3-D spherical surface is covered with an imaginary three-dimensional grid (say 5X5 
X 5) with grid elements being cubes of the same size. For each grid element, an 
associated counter keeps count of the number of field lines created up to that point that 
have intersected that cubical grid element. For each of the grid elements, we create a new 
field line passing through its center only if none of the previously created lines have 
intersected that grid element. The direction of the field along each curved field line is 
indicated through an arrowhead. 

In mode two, the vector field is visualized only on the surface of the selected spherical 
region. The spherical surface is divided into an imaginary, approximately uniform, 
surface grid. The vector field is depicted at each of these points by drawing a vector of 
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length proportional to the strength of the field at that point and in the direetion of the field 


at that point. For the imaginary grid on the spherical surface defined in terms of 0 and cf), 
and for a given r, 0 is incremented by equal amounts, A0, and for a given value of 0, (|) is 
incremented by equal amounts of A0/sin(0). 

Visualization of a Scalar Field 


A scalar field is also visualized in two modes. In mode one, the scalar field is visualized 
through depiction of points with density proportional to the strength of the scalar field at 
the corresponding location in 3-D space. In mode two, small spheres with volume 
proportional to the magnitude of the field are drawn at the centers of the same imaginary 
grid elements as used to control the density of flux lines for the vector field. The 
differentiation between positive and negative values of the field is made through the use 
of two colors one for the positive values of the scalar field and the other for the negative 
values of the scalar field. 

Navigation Through the Fields 

Navigation is accomplished by using the electronic controller/pointer (“wand”); the 
wand as programmed for this application is as follows: 

Wand Button 1: Press and release to change viewing mode 
Button 3: Press and release to reset to initial settings 


8 



Button 2: Press and release to change joystick mode 
joy_stick_mode == 1 : 

joy_stick_x: translates in the direction of the wand 
joy_stick_y: rotates about CAVE's y-axis 
joy_stick_mode == -1: 

Joy_stick_x: To change scale factor 
joy_stick_y: To change the radius of the view sphere 

Visualization of the Earth’s Magnetic Field 

For the computation of the earth’s magnetic field, we use the International Geomagnetic 
Reference Field (IGRF) FORTRAN program to compute the magnetic field at any given 
point in the surrounding 3-D space of the earth. A FORTRAN subroutine can be called 
within a C-program as follows: 

Subroutine in FORTRAN: 

Subroutine subl(x, y, z, arrl, arr2,..) 

Calling the above subroutine in a C-program: 

Subl_(&x, &y, &z, arrl, arr2,. . .), i.e. the name of the subroutine needs to be extended by 
the addition of an (underscore) character, and pointers to call time parameters of the 
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matching type need to be used instead of the parameters themselves. Array names (e.g. 
arrl, arr2) are to be used as they are, because these are pointers to start with. 

In the IGRF subroutine, given an exterior point to earth’s surface in terms of its latitude, 
longitude, and height (lat, long, height), returns the magnetic field at that point as B = 
(BNorth, Be 3 ,j, Bdj,„„). Given the (x, y, z) values in the CAVE’s world coordinate system, 
(lat, long, height) are obtained as follows: 

r = sqrt(x^ + y^ + z^) 

height = r * 2000.0 - 6378.0 Km., where 2000.0 is a scaling factor used 
lat = sin '(z/r) 

long = tan '(y, x), -180°< long < 180° 

if (long < 0.0) long = long + 360, 0° < long < 360° 

The returned value of the magnetic field in the (BNonh, B^ast, Bd„„„) format can be 
transformed to (Bx, By, B^) as follows: 

Bx = - sin(lat) * cos(long) * B^ - sin(long) * B^ - cos(lat) * cos(long) * B^ 

By = - sin(lat) * sin(long) * B^ + cos(long) * B^ - cos(lat) * sin(long) * B^ 

B^ = - cos(lat) * Bn - sin(lat) * B^ 
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Visualization of the Distribution of the Trapped Protons and Electrons in the Earth’s 
Magnetic Field : 

Geomagnetically trapped protons have been specified by the NASA standard APSMIN^ 
database; the electrons use the AESMAX"^ model. These data are consistent with 
maximum observed fiuxes, which occur at different times during the solar cycle, and 
relate directly to the local magnetic field strength and L-shell values. To achieve the 
visualization of the trapped particles, the 3-D space around earth was divided into a three 
dimensional grid of rectangular cells. In each cell, a number of points proportional to the 
inward flux of these particles in that cell was plotted. 

Texture Mapping (of earth’s texture) on a spherical surface from an RGB File: 

The texture on the sphere representing the earth was rendered by mapping an RGB file 
“land_ocean2.rgb” on to the texture file in OpenGL, the graphics modeling and rendering 
software, and then mapping that on the sphere representing the earth through automatic 
generation of corresponding texture coordinates for the surface of the sphere. 

Conclusion 

A system for the visualization of vector and scalar fields in the CAVE has been 
developed. The system has been successfully used for the visualization of the earth’s 
magnetic field, and the visualization of geomagnetically trapped protons and electrons in 
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the space surrounding earth. The future plans for the use and extension of the system 
include the study of the magnetic field around the revolving ISS (International Space 
Station). 
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